Rashba coupling in quantum dots in the presence of magnetic field 
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I. INTRODUCTION 

The Schrodinger equation describing electron in a two-dimensional quantum dot normal to the z axis is of the form 

p2 



2M, 



eff 



Vcix,y) + VR + Vz]^ ^ E^f, (1) 



where Me// is the effective electron mass. The vector potential A = y,x,0) of a magnetic field oriented per- 
pendicular to the plane of the quantum dot leads to the generalized momentum P = p + -A. We have the usual 
expression for the Zeeman interaction 

Vz = \gPBBa,, (2) 

where g represents the effective gyromagnetic factor, //s is the Bohr's magneton The Rashba spin-orbit interaction 
[l|, [2| is represented as 

Vr = aR^a^Py - ayPx). (3) 
The Pauli spin-matrices are defined as standard, 
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A confining potential is usually assumed to be symmetric,Vc(x, — Vc{p),p = \/ + iP- There are two model 
potentials which are widely employed in this area. The first is a harmonic oscillator potential [3, |3| • Such a model 
admits the approximate (not exact) solutions of Eq. (1). The second model is a circular quantum dot with hard walls 
0,0 = for p < Po, = oo for p > Po • This model is exactly solvable. In the framework of above models 

the number of allowed energy levels is infinite for the fixed total angular momentum in the absence of a magnetic 
field. 

In this paper, we propose new model which corresponds to a circular quantum dot with a potential well of finite 
depth: Vc{p) — for p < p^, Vc{p) = V ^ constant for p > po . Our model is exactly solvable and the number of 
admissible energy levels is finite for the fixed total angular momentum in the absence of a magnetic field. The present 
solutions contain, as limiting cases, our previous results [3] (no external magnetic field). 



II. ANALYTICAL SOLUTIONS OF THE SCHRODINGER EQUATION 

The Schrodinger equation (1) is considered in the cylindrical coordinates x = p cos ip,y = psimp. Further it is 
convenient to employ dimensionless quantities 

p 2Me// , _ 2M,ff , 2Afe// ,eBpl g_M^ , . 
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Here Me is the electron mass. As it was shown in [H equation (1) permits the separation of variables 

vl/„(r,^) =M(r)e^™'^ +w(r)e^('"+i)'^ , m = 0, ±1, ±2, . . . (5) 



due to conservation of the total angular momentum + ^(Jz 
We have the following radial equations 



d u 1 du m 2 ^ , 

——r H — h (e — v)u —u — 2bmu ~ b r u — 4sbu 

dr^ r dr 

dw m + 1 

— \ w + brw 

^ dr r 

d'^w 1 dw , , (m+1)^ 2 2 . , 

—TT H f- (e — v)w T. w — Ibyra + \)w — b r w + Asbw 

dr'' r dr 

du m \ 

— — I u + bru . (6) 

dr r J 

In the requirements u{l) = w{l) = were imposed. In our model, we look for the radial wave functions u{r) 

and w(r) regular at the origin r — and decreasing at infinity r —^ oo. 
Following we use the substitutions 



u{r) = exp (^-^^ (x/5r)l"l/(r), wir) = exp (^-^^ (V5r)l'"+i|g(r) 



(7) 



which lead to the confluent hypergeometric equations in the case a = 0. Therefore we attempt to express the desired 
solutions of Eq. (6) via the confluent hypergeometric functions when a ^ 0. 

We consider two regions r < 1 (region 1) and r > 1 (region 2) separately. 

In the region 1 (v — 0), using the known properties 

(/? - 1 - OM{a, P, + S.^^^^^^ = {P- l)M{a - 1, /3 - 1, (8) 

of the confluent hypergeometric functions M(a, (3, ^) of the first kind ^8|| it is easily to show that the suitable particular 
solutions of the radial equations are 



where 



for m = 0,1,2. 



uiir) = exp (^-^^ (V5r)l"l (ci_/i_(r) + ci+/i+(r)) , 
w,ir) = exp (^^) (V6r)l™+i| (^^) (ci^gi-(r) + ci+gi+(r)) , (9) 



fi^{r) = M(m + 1 - fcjF,m + l,6r^), 
,' kf \ M(m + 1 - kf,m + 2,br^) , , 

^^-^^^ ^ ( (^j (4^+(4.)^^:-H.-i/2) 



/i^(r) =Af(l-fcf,_,„ + l^for2)^ 

. , M(-kf,-m,br'^) , , 

gi^r =m — ^ ^— '- (11 

■^^^ ^ (-fcf + (46)-ie + .s-l/2) ^ ^ 



for m = —1, —2, —3... and 



1 / a' I a? f^b^ ^ 



^i^ = ¥h+y^"\/'+T + (-) (--V2)M- (12) 
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Here ci_ and ci+ are arbitrary coefRcients. The functions ui{r) and z«i(r) have the desirable behavior at the origin. 
In the region 2 (v > 0), using the known properties 



Uia, (3, - = Uia, 13+1,0, 



-U{a-l,P-l,0 



(13) 



of the confluent hypergeometric functions U{a, (3,0 of the second kind it is simply to get the suitable real solutions 
of the radial equations: 



'W2{r) ~ exp 



U2{r) = exp 



(\/6r)'"l (c2-/2-(r)+C2+/2+(r)) 



where 

for m = 0,1, 2..., 



([/(m + 1 - fcj , m + 1, 6r^) =F C^(™ + 1 - ,m + 1, fer^)) , 



;7(m + 1 - ,m + 2, for^) 



U{m + 1~ ,m + 2,br'^) 



{-k^ +{Ab)-^{e-v) + 8-1/2) (-fc2+ + (46)-i(e-w) + s-l/2) 



{U{1 ~ k2,-m + 1, br^) T U{1 - kt , -m + 1, br'^)) , 



U{—k2 , — m, br^) 



T 



U{—k^, —TO, br^) 



i-k^ + {Ab)-^{e-v) + s-l/2) i-k+ + {Ab)-^{e-v) + s-l/2) 



for 



-1, —2, —3... and 



T 



(.-1/2)^ 



(14) 



(15) 



(16) 



(17) 



Here C2_ and C2+ are arbitrary coefficients. The functions U2{r) and W2{r) have the appropriate behavior at infinity. 
We assume the realization of condition 



e < V* = V 



46 



(,s-l/2)^ 



(18) 



which means that electron belongs to a quantum dot. We can also obtain the exact solutions when e > w*. However, 
in this case we cannot consider electron as belonging to a quantum dot. 
The continuity conditions 

ui(l)-U2(l) = 0, u;i(l)-u;2(l) = 0, u'i(l) - u^(l) = 0, w'l (1) - ^(1) = (19) 

for the radial wave functions and their derivatives at the boundary point r — 1 lead to the algebraic equations 



TA(m, e, v.a, b, s 



ci+ 

\ C2+ 







(20) 



for coefficients ci_,ci+,C2- and C2+ where 
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(21) 
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Hence, the exact equation for energy e(m, v, a, b, s) is 

detT4{m,e,v,a,b,s) =0. (22) 



This equation is solved numerically. 
The desired coefficients are 



ci+ \ / -/i-(l) 



C2- 



Ci-T^\m,e,v,a,b,s)\ -9i-{l) , (23) 



where 



C2+ / V -/^(l) 



//l+(l) -/2-(l) -/2+(l) 

Ts{m,e,v,a,b,s)= gi+{l) -52-(l) -52+(l) 1 • (24) 

V/(+(i) 

The value of ci_ is determined by the following normalization condition Jg°° (u^(r) + w^(r)) rrfr = 1. 

III. NUMERICAL AND GRAPHIC ILLUSTRATIONS 

Now we present some numerical and graphic illustrations in addition to the analytical results for the ground and 
first excited states in the particular cases to = 1, m = — 2 at fixed s = 0.05. 

Tables show the energies e for diS^erent values of the Rashba parameter a, the well depth v and the magnetic field 

b . 

Figures demonstrate the examples of continuous radial wave fmictions for v = 100, a — 2,b = 5. Solid lines 
correspond to the functions u{r) and dashed lines correspond to the functions u'(r). We see that the radial wave 
functions rapidly decrease outside the well. The values of coefficients are ci_ = —0.443198, ci+ = 5.63821, C2- = 
-20148.7 and C2+ = 6147.72 in the case of Fig. 1 and ci_ = -0.065795, ci+ = -2.5723, C2- = -31387 and C2+ = 
2251.07 in the case of Fig. 2. 



TABLE I: Energy levels for m = 1. 



v = 50 


V = 100 


b 


a = 1 




a = 2 


a = 1 


a = 2 







10.23 


20.31 


7.97 20.73 


11.16 22.08 


8.85 


22.59 


0.5 


11.33 


22.50 


8.84 23.09 


12.25 24.24 


9.73 


24.94 


1 


12.65 


24.95 


9.92 25.70 


13.55 26.64 


10.81 


27.53 


1.5 


14.18 


27.68 


11.23 28.52 


15.05 29.29 


12.08 


30.32 


2 


15.93 30.67 


12.74 31.54 


16.74 32.18 


13.56 


33.32 


2.5 


17.88 




14.47 34.68 


18.62 35.30 


15.23 


36.47 


3 






16.39 37.92 


20.68 38.65 


17.08 


39.75 


3.5 






18.49 


22.90 42.20 


19.10 


43.11 


4 






20.75 


25.27 45.94 


21.28 


46.47 


4.5 






23.16 


27.77 


23.61 


49.82 


5 






25.70 




26.07 


53.17 


5.5 










28.65 


56.57 


6 










31.32 


60.06 



5 



TABLE II: Energy levels for m = -2. 



v^50 


V = 100 


b 


(I = 1 




a = 2 




(/ = 1 




a = 2 







10.23 


20.31 


7.97 


20.73 


11.16 


22.08 


8.85 


22.59 


0.5 


9.36 


18.41 


7.33 


18.63 


10.26 


20.18 


8.18 


20.47 


1 


8.71 


16.79 


6.89 


16.79 


9.57 18.52 


7.70 


18.60 


1.5 


8.26 


15.45 


6.66 


15.21 


9.08 


17.12 


7.40 


16.98 


2 


8.02 


14.37 


6.61 


13.91 


8.77 


15.96 


7.29 


15.59 


2.5 


7.96 


13.56 


6.74 


12.87 


8.63 


15.05 


7.33 


14.45 


3 


8.08 


12.99 


7.01 


12.08 


8.66 


14.36 


7.53 


13.55 


3.5 


8.34 




7.42 


11.53 


8.84 


13.90 


7.86 


12.87 


4 






7.94 


11.21 


9.15 


13.64 


8.31 


12.41 


4.5 






8.55 


11.09 


9.58 


13.58 


8.86 


12.15 


5 






9.24 


11.18 


10.11 


13.70 


9.50 


12.08 


5.5 






9.98 


11.43 






10.20 


12.19 


6 






10.77 


11.84 






10.95 


12.46 





FIG. 2: Radial wave functions for m = 



-2,e = 12.0827. 
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IV. CONCLUSION 



So, we have constructed new exactly solvable and physically adequate model to describe the behavior of electron 
in a semiconductor quantum dot with account of the Rashba spin-orbit interaction and the external magnetic field. 
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